Draft version February 2, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



O 
O 
(N 

> 

O 

(N 



I 

o 



en 

> 

in 

00 

o 
m 



o 



X 



A ROTATING HOLLOW CONE ANISOTROPY OF TEV EMISSION FROM BINARY SYSTEMS 

A.Neronov 

INTEGRAL Science Data Center, 16 ch. d'Ecogia, CH-1290, Versoix, Switzerland 

AND 

M . C HERN YAKOVA 

Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, Dublin 2, Ireland 

Draft version February 2, 2008 

ABSTRACT 

We show that TeV 7-ray emission produced via interactions of high-energy particles with anisotropic 
radiation field of a massive star in binary systems should have a characteristic rotating hollow cone 
anisotropy pattern. The hollow cone, whose axis is directed away from the massive star, rotates with 
the period equal to the orbital period of the system. We note that the two maxima pattern of the TeV 
energy band lightcurve of the 7-ray loud binary LS 5039 can be interpreted in terms of this rotating 
hollow cone model. Adopting such an interpretation, we are able to constrain the geometry of the 
system - either the inclination angle of the binary orbit, or the elevation of the 7-ray emission region 
above the orbital plane. 
Subject headings: gamma rays: theory — radiation mechanisms: non-thermal — binaries: general 



Introduction. 7-ray-loud binary systems are a newly 
identified class of sources in which either accretion onto 
the compact object (a neutron star, or a black hole), 
or interaction of an outflow from the compact object 
with the wind and radiation from a massive companion 
star leads to the production of very-high energy (VHE) 
7-ray emission. Three such systems, PSR B1259-63, 
LS 5039 and LSI -1-61 303, have been firmly detected 
as persistent or regularly va r iable TeV 7-ray e mitters 
(|Aharonian et al.ll2005l [200I [Albert et al.ll2006D . The 
VHE 7-ray emission from the 7-ray-loud binaries is vari- 
able on the orbital period (or shorter) time scale. This 
implies that the emission region is located close to the bi- 
nary system, in a highly inhomogeneous and anisotropic 
particle and photon background produced by massive 
companion star. 

In what follows we show that if the 7-ray emission from 
such a region is produced in interactions of isotropically 
distributed VHE particles with photons from the mas- 
sive star, it should have a characteristic "rotating hollow 
cone" anisotropy, i.e. most of the photons are emitted 
at a certain angle Co with respect to a symmetry axis 
directed radially away from the massive star. Orbital 
motion of the emission region around the massive star 
leads to the rotation of the emission cone. Rotation of the 
hollow cone on the orbital time scale leads to the appear- 
ance of 0, 1, or 2 maxima in the orbit-folded lightcurve 
in the VHE band, occurring at the phases when the line 
of sight is inclined at an angle Co with respect to the 
cone axis, i.e. at the moments of passage of the of the 
cone through the line of sight (similarly to the hollow 
cone models of period-folded light curves of pulsars, see 
e.g. (jLvne & Craham-Smithl[2005[ )V 

The orbital modulation of the 7-ray flux, related to 
the passage of the hollow cone through the line of sight 
could be most clearly detected if there are no additional 
sources of the modulation, related e.g. to the ellipticity 
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of the binary orbit, absence of spherical symmetry of the 
wind/radiation from the companion star etc. Among the 
three 7-ray-loud binary systems mentioned above, the 
system LS 5039 is characterized by the lowest ellipticity 
of the orbit. In this system the compact object orbits a 
06.5V star which emits isotropic stellar wind (contrary 
to the other two systems in which the massive star is of 
the Be type). 

The influence of the anisotropy of the photon field of 
the massive star on the properties of the 7-ray emission 
in binaries in general, and in LS 50 3 9 in particu lar, was 
first studied by Khang ulvan et all ()2005i |2007| ). Here 
we calculate the angular brightness profile of the hollow 
cone in LS 5039, and find that the observation of the two 
maxima of the orbit-folded lightcurve constrains the in- 
clination of the binary orbit to be i > 40°, if the emission 
is produced in the vicinity of the compact object. This 
result can be stated also in an opposite way: if the in- 
clination of the binary orbit is i < 40° , the two maxima 
structure of the orbit-folded lightcurve can be explained 
only if the VHE 7-ray emission region is displaced from 
the position of the compact object. This can be the case 
if the emission is produced in a jet. In this latter case, 
we show that the existence of the two-maxima of the 
lightcurve constrains the elevation of the emission point 
above the orbital plane. 

Anisotropy of VHE ^-ray emission in a central photon 
field. Consider the 7-ray emission produced by inter- 
actions of VHE particles X (e.g. protons or electrons) 
with the soft photon field in the vicinity of a massive 
star. Assume for simplicity that the size of the emission 
region is much less than the distance from the region to 
the center of the star and that the VHE particles in the 
region have isotropic velocity distribution. In spite of 
the isotropy of the VHE particle distribution, the 7-ray 
emission will be anisotropic. The anisotropy arises be- 
cause of the Doppler effect which leads to the decrease 
(increase) of the rate of interaction of the VHE particles 
co-moving with (moving oppositely to) the soft photon 
field of the massive star. 
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The interaction rate of particles X with monienta 
Pjc with soft photons w ith momenta p* is given by 
(iLandau fc Lifshitzlfl980l) 
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where a is the interaction cross-section, nx(Px) is the 
particle distribution, n*(p*) is the soft photon distribu- 
tion and Px ■ P* is the scalar product of the 4-momenta 
of the interacting particles, 

Px-p*=^^.e*(l-cosC) , (2) 

(C is the angle between the particle velocities) which con- 
tains the Doppler factor, (1 — (3cos(^) (we assume that 
particle velocity is /3 ~ 1). 

If the soft photon field in the emission region would be 
isotropic, integration of the interaction rate of Eq. ([T]) 
over the soft photon angular distribution would average 
out of the angular factor (1 — cos^), so that the 7-ray 
emission intensity, which is proportional to the interac- 
tion rate, would be isotropic. However, since all the soft 
photons crossing the emission region move in the same 
direction (away from the massive star), the interaction 
rate depends on the angle ^ between the direction of 
emission and the direction " from the massive star" . Par- 
ticles X comoving with the soft photon field (moving at 
the angles C ^ 0°) interact more rarely than the particles 
moving opposite to the photon field (at C -^ 180°). 

The reduction of the interaction rate has two- fold con- 
sequences. On one hand, it leads to the reduction of the 
power of the 7-ray emission by the particles moving in 
the direction away from the star. On the other hand, the 
reduction of the interaction rate of the soft photons with 
the emitted VHE 7-rays facilitates the escape of 7-rays 
moving in the direction away from the star. A competi- 
tion between the decrease of the 7-ray production rate, 
7?-prod(C)i and the increase of the 7-ray "survival prob- 
ability" (i.e. of the exp(— Tabs(C))i where rabs(C) is the 
optical depth with respect to the pair production) leads 
to the appearance of a maximum of the 7-ray flux 

i^^(C)~7^p,od(C)e"""^=^'^^ (3) 

at an angle < Co < tt, i.e. to the appearance of a hollow 
cone anisotropy pattern. 

At large distances from the massive star one can ap- 
proximate the angular distribution of the soft photons by 
that of a point source. Under this simplifying assump- 
tion, one finds that the integration over the angular dis- 
tribution of the UV photons is easily performed and the 
resulting expression for the rate of production of 7-rays 
takes the form 

''^prod(C) - "-Xn*tTprod(l " COS C) , (4) 

where CTprod is the cross-section of production of 7-rays in 
interaction of the particles X with the soft photon field. 
An estimate of the optical depth for the 7-rays, es- 
caping from the production region, can be obtained by 
multiplying the absorption rate per 7-ray on the size of 
the absorbing region. The absorption rate is given by 
the expression ([1]) in which the particle X is a 7-ray and 
the cross-section a is the pair production cross section, 
(Tabs- Estimating the size of the absorption region to be 
of the order of the distance D of the emission point from 
the massive star, one finds 

Tabs ^ TZabsD/rij ~ n*crabs£>(l - cosC) . (5) 




Fig. 1. — Thick solid curve: angular brightness profile, Eq. 
(|6]l, in the approximation of constant cross-sections. Thin dotted 
and dashed curves show the 7-ray production rate and absorption 
coefficients as functions of the angle (J. 




Fig. 2. — Model of the anisotropy of the 7-ray emission. Maxima 
of the lightcurve are expected during the passage of the line of 
sight through the hollow cone, whose axis is directed along the line 
connecting the emission region. 



Substituting Eqs. (g]),® into Eq. Q one finds 



-^7(0 ^ crprod(l - cosC)e 



-n,£l(Tabs(l-cosC) 



(6) 



Two effects affect the anisotropy pattern of the 7-ray 
emission. First, the explicit dependence of F-y on cos^ 
is introduced in Eq. ^ by the Doppler effect. An ad- 
ditional implicit dependence on ^ is introduced through 
the energy dependence of the interaction cross-sections 
Cprodj^abs- Indeed, in general a — (t{Ecu), where the 
Eqm if the center-of-mass energy, which, in the case 
^CM ^ 'mx depends on C as Eqm ~ -y/l ~ cos^. 

The anisotropy pattern resulting from the Doppler ef- 
fect can be found if one ignores the energy dependence 
of CTprod, cTabs- This is donc in Fig. [1] (the constant 
To — n^ZJcTabs is taken to be tq = 3). From this figure 
one can see that most of the 7-ray flux is emitted along 
a "thick hollow cone" with the opening angle (0 and the 
thickness comparable to the opening angle, A^ ^ Co- 
The energy dependence of the interaction cross-sections 
leads to the energy-dependent angular brightness profile 
of the thick hollow cone. The cone becomes wider at 
higher energies, where the absorption is less efficient, see 
Fig. [5] below. The non-zero brightness at C = is due to 
the flnite radius of the star. 

Geometrical model of variability ofj-ray emission. If 
the location of the emission region is determined by the 
position of the compact object (e.g. the 7-ray emission is 
produced in the vicinity of the compact object, or in the 
jet emitted by the compact object), the orientation of 
the hollow cone changes when the compact object moves 
around the star, as it is shown in Fig. [2] (where the 
emission region is supposed to be situated at an elevation 
h above the compact object, so that the hollow cone axis 
is inclined at an angle x — atan(/i/£>), where D is the 
binary separation distance). 

Changes of the orientation of the hollow cone with re- 
spect to the line of sight should lead to the orbital mod- 
ulation of the observed 7-ray flux. Depending on the re- 
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Fig. 3. — A phenomenological model fit to the Ughtcurve of LS 
5039. The model consists of a constant plus two gaussians of equal 
width at equal distance from the inferior conjunction (shown by a 
vertical line). 

lation between the inclination angle i of the binary orbit 
plane, and the opening angle of the hollow cone, ^Oj two 
characteristic patterns of the orbital modulation are pos- 
sible. If the inclination of the orbit is i > 7r/2 — Cq — Xj the 
direction of the line of sight passes through the "wall" of 
the cone two times per orbit. This should lead to the ap- 
pearance of two maxima in the orbit-folded lightcurve. If 
the binary orbit is circular, the maxima of the lightcurve 
occur when the true anomaly of the orbit (angular or- 
bital phase < $ < 360°, counted from the focal point 
of the orbit, $ = 0° at the periastron) takes the values 

$1,2 = ^-inf ± A$ (7) 

where $inf is the anomaly of the inferior conjunction and 

A<i> = acos [1 — (sin(i -I- x) ~ cosCo)/(cosxsini)] (8) 

If the binary orbit is elliptical, an additional orbital mod- 
ulation of the 7-ray flux can occur because of the varia- 
tion of the distance of the emission region from the mas- 
sive star (which leads to the modulation of the 7-ray 
production/absorption rates) with the orbital phase. As 
the inclination of the orbit decreases to i < 7r/2 — Co ~ Xj 
the two maxima at $1 and <i>2 merge at the phase $inf- 

Since the opening angle of the hollow cone, Co depends 
on the 7-ray energy, the condition i > 7r/2 — Co (£'7) — X 
can be satisfied only in a certain interval of energies, so 
that a merger of the two maxima <I>i^2 = 'I'inf ± A$(_E-y) 
at the phase of inferior conjunction can be observable at 
a particular energy Eq at which Co(-£'o) = 7''/2 ~ i — x- 

The case of LS 5039. LS 5039 is one of the sev- 
eral X-ray bi naries detected a.s sou rces of the VHE 7- 
ray emission (jAharoni an et al.]l2006[ ). In this binary sys- 
tem the compact object rotates with the period P — 
3.9078 ± 0.0015 d around a massive 06.5V star. The or- 
bit is eccentric (e = 0.35). The inc lination of the orbit is 
20° < i < 60° (|Casareset al.l[2005l) . 

The orbit-folded lightcurve of the source at the energies 
E > 1 TeV (Aharonian ct al. 2006) has two maxima at 
the orbital phases (pi ~ 0.55, 02 — 0.85. The two maxima 
are apparently not symmetric: the first maximum around 
(j)i spans a broader range of the orbital phase, while the 
second maximum is more sharp. The second maximum 
happens closer to the phase of the inferior conjunction, 
(/)inf — 0.716, than the first one. Replotting the orbit- 
folded TeV 7-ray lightcurve as a function of the true 
anomaly, $, (to produce Fig. [H we have calculated the 
true anomalies of each data point of the top panel of Fig. 
5 of I Aharonian et all (|2006) and rcbinncd the data into 
bins of the width d$ — 15°) rather than as a function of 
(p , a symmetry in the positions of the two maxima can 
be found. 
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Fig. 4. — Dependence of the opening angle of the hollow cone on 
the inclination of the binary orbit in LS 5039. Grey shaded region 
indicates the range of C,o found from numerical calculations. 

Namely, the lightcurve can be satisfactory fitted with 
a phenomenological model which is a sum of a constant 
(weakly modulated emission which can be produced e.g. 
at larger distances) plus two gaussians, whose centers are 
equally spaced from the phase of the inferior conjunction, 
$inf ~ 224° (see Fig. ^. The positions of the centers 
of the gaussians, found by the fit, are $i_2 — ^inf ± 
35°, while the widths of the gaussians are nearly equal, 
(5$i,2 — 22.5° (if the phase of the inferior conjunction 
is left free, while fitting, the fit finds the phase $inf — 
226.5°, c onsistent with the va lue $i„f ~ 224.2° ± 3.3° 
found bv ICasares et al.l ()2005[ )). The gaussian centered 
at the phase $2 is found to have ~ 1.5 times higher 
normalization than the gaussian centered at ^i. 

The observed symmetry of the positions of the max- 
ima of the orbit folded lightcurve can be interpreted, in 
a straightforward way, in terms of the " rotating hollow 
cone" model, discussed in the previous sections. In par- 
ticular the phase shift of Eq. © is A$ ~ 35°. From Eq. 
([5]) (which can be used, for a low-eccentricity orbit, as a 
first approximation) one can find a relation between i and 
Co, shown in Fig. 31 Taking into account the constraint 
on th e inclination angle, 20 < i < 60 (jCasares et al.l 
l2005f ). we find that the opening angle of the hollow cone 
anisotropy pattern is constrained to 40° < Co < 75°, if 
the emission is assumed to come from the vicinity of the 
compact object (x — 0° in Fig. [2). 

A constraint on the geometry of LS 5039. The opening 
angle of the hollow cone, Co (£-7)7 can be found, once the 
emission process leading to the 7-ray production and the 
location of the emission region are known. Comparing 
the theoretically predicted value of Co (^7) to the one 
implied by the data one can, in principle, constrain the 
geometry of the system. In particular, assuming that 
the location of the emission region is known, one can 
constrain the inclination of the binary orbit with respect 
to the line of sight. Otherwise, if the inclination of the 
orbit would be known, one would be able to constrain the 
location of the emission region, in particular, its distance 
from the star and/or elevation above the orbital plane. 

Different locations of the VHE 7-ray emission region 
are assumed in different models of activity of LS 5039. 
In the model of "co mpact pulsar wind nebula" (see e.g. 
(jDubus et al.ll2007[) ) the emission region is assumed to 
surround the compact object (/i = in notations of Fig. 
[2). In a " microquasar" model the TeV emission is as- 
sumed to be produced in a jet, so that the emission re- 
gion is displaced from the position of the compact object 
(it is situated above or below the orbital plane, h^ 0, see 
e.g. ([Khangulvan et aT] |2007')V In both types of models 
the VHE 7-rays are produced via the inverse Compton 
scattering of the soft photons by the VHE electrons. 
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Fig. 5. — Angular brightness profile of the hollow cone in LS 
5039, found assuming that emission is produced at the distance of 
periastron (black) or apastron (grey) of the orbit. Thin dotted and 
dashed curves show the production rate and absorption coefficient 
for an emission region at the periastron distance. 

To find the rate of production of 7-rays via the inverse 
Compton scattering, one has to numerically integrate the 
Eq. H]), with CTprod being the Klein-Nishina cross section, 
over the angular and energy distributions of soft photons 
at a distance D from the massive star. The result of 
such integration is shown by the dotted curves in Fig. 
[5] for electron energies. Eg = 1 and 10 TeV for the case 
when the emission region is situated at the distance of 
the periastron of the binary orbit. 

To find the optical depth Tabs (C) for the 7-rays emitted 
at different angles C (see Eq. [3]), one first calculates the 
7-ray absorption rate, given by Eq. ([T]) with the cross 
section a being the pair production cross-section, at each 
point of the 7-ray trajectory. This is done similarly to 
the calculation of the inverse Compton scattering rate, 
via an integration over the soft photon angular distri- 
bution. Next, one has to integrate the absorption rate 
along the 7-ray trajectory, from the emission point to in- 
finity. The result of such numerical integration is shown 
by the dashed curves in Fig. [5] for 7-ray energies E^ — 1 
and 10 TeV, assuming that emission is produced at the 
distance of the periastron of the binary orbit. 

Substituting the numerically found 7?.prod(C7 ^7) and 
■'abs(C: ^7) into Eq. ^ one can find the angular bright- 
ness profile of the hollow cone for different orbital phases 
and different photon energies. The angular brightness 
profile, calculated for the periastron/apastron of the bi- 
nary orbit, is shown by the thick, solid, black/grey curve 
in Fig. [5] We assume that the energies of 7-rays are 
approximately equal to the energies of the primary elec- 
trons, which is true in the Klein-Nishina regime of the 
inverse Compton scattering. 

The first maximum of the TeV band lightcurve of LS 
5039 takes place at the orbital phase $1 ~ 224° - 35° ~ 



193°, close to the apastron of the orbit. From Fig. [5] 
one can see that at this phase the opening angle of the 
cone changes in the range 30° < (o{E^) < 60° when the 
7-ray energy changes from 1 to 10 TeV. Comparing this 
numerically calculated range of values (shaded region in 
Fig. HI to the one implied by the observational data, one 
can find from Fig. 21 that if the VHE 7-ray emission is 
produced close to the compact object (i.e. at x = 0°) the 
inclination angle of the binary orbit should be i > 40°. 

The constraint on i,x can be re- formulated in a differ- 
ent way: if th e inclination a n gle of the orbit is small, e.g. 
i ~ 25° (see ICasares et al.l ()2005[ )). the 7-ray emission 
is not produced close to the compact object. Instead, 
from Fig. [Hone can find that in this case the elevation 
X of the emission region above the orbital plane should 
be X > 15°. If the emission region is located in a jet-like 
outflow orthogonal to the orbital plane, the emission re- 
gion should be situated at the height h = D tanx > 0.3D 
above the orbital plane, where D is the binary separation 
distance (see Fig. [2] for notations). 

Summary. We have shown that VHE 7-ray emis- 
sion from 7-ray-loud binaries is expected to have a ro- 
tating hollow cone anisotropy pattern, determined by 
the Doppler effect in the anisotropic radiation field of 
a massive star (Fig. [Ij . This anisotropy leads to the ap- 
pearance of a double-peak structure of the orbit-folded 
lightcurve, with the two peaks situated at equal distance 
A<I> (Eq. [5]) from the phase of the inferior conjunction. 
We have demonstrated that such a symmetric double- 
peak structure is observed in the particular case of LS 
5039. In this case, a measurement of the phase shift A<I> 
enables to find a relation between the opening angle of 
the hollow cone, Co, and the inclination of the binary or- 
bit, i (see Fig. E]). Comparing the value of Co, inferred 
from the data, to the one found from numerical calcula- 
tion of the angular brightness profile of the cone (Fig. \5^ 
we were able to constrain the inclination of the binary 
orbit and/or the elevation of the VHE 7-ray emission 
region above the orbital plane in this particular source. 

In the particular case of LS 5039, the rotating hollow 
cone model discussed above can be tested if the statistics 
of the signal from the source becomes high enough to 
allow a splitting of the source lightcurve at the energies 
E > 1 TeV onto two energy bins (e.g. 1 TeV< E < 
10 TeV and £■ > 10 TeV). In this case the predicted 
shifts of the two maxima of the lightcurve toward each 
other (or even a merger of them) at the higher energies 
should be observable. If observed, such an effect would 
be a clear evidence in favor of the proposed model. 

We thank F.Aharonian and A.Zdziarsky for the discus- 
sions of the subject. 
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